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■ Results of high-resolution simultaneous multi-frequency 8.1-15.4 GHz VLBA polarimetric ob- 

\^ i servations of relativistic jets in active galactic nuclei (the MOJAVE-2 project) are analyzed. We 

compare characteristics of VLBI features with jet model predictions and test if adiabatic expan- 
sion is a dominating mechanism for the evolution of relativistic shocks in parsec-scale AGN jets. 



OO ' We also discuss magnetic field configuration, both predicted by the model and deduced from 

o 



electric vector position angle measurements. 
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1. Introduction 

In spite of the fact that our understanding of active galactic nuclei phenomenon has been 



significantly improved over the past decade [11], there are many aspects of nature of the relati vis- 
tic outflows (jets), like their composition, formation, acceleration, coUimation, etc, which remain 
to be strongly debated. VLBI observations provide us with the unique tool to explore the jets in 
AGN with milliarcsecond angular resolution corresponding to parsec-scale resolution in linear size. 
Typically, parsec-scale flows are characterized by one-sided knotty structure with sometimes pro- 
nounced curvature, rapid variations of flux density (e.g. U), and detection of superluminal 
motions [^. 

Relativistic shocks propagating down the jets are expected to be prominent on these scales. 



which is confirmed by detection of strong polarization [ |12| , |13| ] and enhancement of the magnetic 
field in the brightest jet features (from the turnover frequency images) |^. The general evolution 
of a shock-induced flare in a jet is considered by ||^ to pass through the Compton, synchrotron, and 
adiabatic stages, referring to the dominant energy loss mechanism. In the current paper we present 
recent results of testing adiabatic expansion of relativistic shocks in several active galactic nuclei. 

Throughout the paper the ACDM cosmological model with//o = 70 km s^' Mpc^\ H,,, = 0.3, 
and Q.\ = 0.7 is adopted. Spectral index a is defined as Soc v". 



2. Observational data in use 

We use experimental data collected and processed in the framework of the MOJAVE-2^ project. 
The observations were performed by the MOJAVE team in 2006 at VLBA, a VLBI network of ten 
25-m telescopes, in a dual-polarization mode, quasi-simultaneously at four frequencies of 8.1, 8.4, 
12.1, & 15.4 GHz. The scan lengths were chosen to achieve roughly the same image rms at each 
observing frequency. The observed sample consists of 192 sources, see for details [||]. 



3. Results and Discussion 

Using the shock-in-jet-model |Q] we have investigated an evolutionary scenario in which a 
bright jet feature is treated as an independent, geometrically thin relativistic shock developing 
in a flow with the emission dominated by adiabatic energy losses. It is also assumed that the 
jet plasma has a power-law particle energy distribution N{E)dE oc E^'' dE. Following [|l[ and 
applying Gauss's flux theorem and Lorentz transformations to the cross section of the jet, the 
magnetic field can be approximated as B oc J^", where d is the transverse jet size and n describes 
the orientation of the magnetic field (n = I for the transverse field and « = 2 for the longitudinal 
field). Assuming non-changing Doppler factor along the jet, the model brightness temperature 



of each jet component, Tb jet, can be related, as it was shown in [IC], to the measured brightness 
temperature of the core, Tb core, as Tb jet = 7b core('^jet/'^core)^^, where d represents the measured 
sizes of the core and jet components, and ^ =n + l — a(?i + 4/3). 



'Monitoring Of Jets in Active galactic nuclei with VLBA Experiments, 

http : //www . physics . pur due . edu/MO JAVE/ 
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We compare the measured brightness temperatures with those predicted by the model, deter- 
mining the spectral index of each jet component as well as the B-field orientation (if polarization 
emission is detected, preferably at 15.4 GHz, by rotating electric vectors by 90° due to the optically 
thin synchrotron radiation detected in the jets). 

In this paper we discuss the results for two sources, 1128—047 and 2155—152, having promi- 
nent outflows with sufficient number of jet components suitable for the analysis. 



3.1 1128-047 



This flat-spectrum radio galaxy with redshift z = 0.266 (3.91 pc/mas) and mass of the central 
object of 5.2 x 10^ M© rtTSl] shows one-sided jet structure (Fig. left) which was fitted by seven 



distinct Gaussian components at 8.1 GHz. The model fitting was performed in DIFMAP [15]. In 
case of applying an elliptical Gaussian component for the VLBI core the minor axis of the ellipse 
has been used to estimate d, since the VLBI core component has a strong tendency to be elongated 
in the innermost jet direction [Q]. 

Using the VLBI images at 8.1 and 12.1 GHz, we constructed spectral index distribution and 
its sHce along the jet ridge fine (Fig. [l], right). Since the linear polarization in this source has not 
been detected down to about 1 mJy level, we calculated the model brightness temperatures for 
aU of the jet components assuming (i) parallel (Fig. |2[ left) and (ii) perpendicular magnetic filed 
configuration (Fig. ^ right). In the latter case, the values of r^^*"* and rj^°'^ agree within the errors 
that allows us to suggest predominantly transverse orientation of magnetic field in the jet. 




10 5 -5 10 5 
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Figure 1: Left panel: Total intensity image of radio galaxy 1128—047 at 8.1 GHz at epoch of 2006.92 
with the model of 7 Gaussian components superimposed. The lowest contour is plotted at 0.35% of the 
peak brightness of 501 mJy/beam. Shaded ellipse represents the FWHM of the restoring beams of 2.54 x 
1 .04 mas at PA = —4.0°. Right panel: Spectral index distribution in 1 128—047 calculated between 8.1 and 
12.1 GHz with the 12.1 GHz total intensity contours overlayed. White dots represent the total intensity ridge 
line along which we plot the profiles of total intensity and spectral index as an inset image. The lowest 
contour is plotted at 0.84% of the peak brightness of 531 mJy/beam. Shaded ellipse represents the FWHM 
of the restoring beam of 2.06 x 0.89 mas at PA = - 1 .4°. 
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Figure 2: Left panel: Brightness temperature changes along the jet in 1128—047 at 8.1 GHz. Filled blue 
circles are the measured values. Open red circles are the predicted brightness temperatures from the shock- 
in-jet model with adiabatic losses dominating the radio emission assuming predominantly longitudinal mag- 
netic field orientation to the local jet direction. The errors are given at 3(7 level. Right panel: The same 
data as on the left but assuming transverse magnetic field orientation. Big discrepancies (left) and remark- 
able agreement (right) between measured and model brightness temperature values lead to the conclusion of 
perpendicular orientation of the magnetic field in the jet. 



3.2 2155-152 

This flat-spectrum source with redshift z = 0.672 (7.03 pc / mas) and mass of the central object 



of 3.9x lO'^M© [fUl was classified as a BL Lac object by Craine et al. and then as a quasar by 
Veron-Cetty & Veron JI^ on the basis of its optical spectrum. 

Using the parsec-scale spectral index information (Fig. ^ right) taken from the I-maps at 8.1 
and 15.4 GHz (Fig. ^, left) together with the results from model fitting, and also assuming magnetic 
field to be perpendicular to the jet, we have calculated the values of T^"^ and The brightness 
temperatures predicted by the model agree well with those from the data for all the components 
(Fig. ^, right), except the optically thin J2 feature {aj2 = —0.75) located at 3.36 mas from the core, 
for which = jobs^jmod 20. There are two possible ways to explain this discrepancy: 

(i) The J2 component is the region of the outflow where the particles are reaccelerated through 
the interaction with the surrounding medium; an evidence for this is the formation of a shear layer 
as is described in the decelerating jet model [^: in the central channel, where magnetic field is 
predominantly tranverse, plasma moves with relativistic speed but slows down at the edges due to 
interaction with an ambient medium making a sheath of longitudinal magnetic field registered on 
polarization images (Fig. ^ left). 

(ii) The Doppler factor varies in the J2 region due to changes of the speed or orientation of the 
jet. The apparent speed jSapp j5 = (10.8 ± 1.3)c is determined from the MOJAVE kinematics |8|]. 
The J2 component is quasi-stationary (during a period of 5.2 yr centroid position was wiggling 
within 0.3 mas in non-radial direction). Then assuming the viewing angle Omax ~ 5.3° (0 ~ l/j8) 
to be constant and applying analytical adiabatic approximation [|I|] for a relativistic jet with pre- 
dominantly transverse B-field, we determined the change of speed Ap < 0.002 between J5 and J2 
regions. Therefore, the detected excess of the Doppler factor 5j2/8j5 = «^i/(3-aj2) = 2.22 can be 
explained by changing the viewing angle that becomes less than the opening angle of the jet ^ 1°. 
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Figure 3: Left panel: Total intensity image of quasar 2155 — 152 at 15.4 GHz at epoch 2006.92 with model 
fitted Gaussian components superimposed. The lowest contour is plotted at 0.24% of the peak brightness 
of 796 mJy/beam. Shaded ellipse represents the FWHM of the restoring beam of 1 .32 x 0.50 mas at PA = 
—5.0°. Right panel: Spectral index distribution in 2155 — 152 calculated between 8.1 and 15.4 GHz with the 
15.4 GHz total intensity contours overlayed. White dots represent the total intensity ridge line along which 
we plot the profiles of total intensity and spectral index as an inset image. The lowest contour is plotted 
at 0.45% of the peak brightness of 1000 mJy/beam. Shaded ellipse represents the FWHM of the restoring 
beam of 1.99 x 0.76 mas at PA = -3.2°. 
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Figure 4: Left panel: Linear polarization images of 2155 — 152 at 15.4 GHz with electric vectors superim- 
posed on epochs 2006.92 and 2008.70. The bottom P-contour is plotted at 1 .8 mJy for both epochs. The blue 
contour represents the total intensity bottom contour at a level of 1 .9 mJy for both epochs. The sheath-spine 
structure with transverse B-field in the central channel of the jet and longitudinal B-field at the edges is seen 
on epoch 2008.70. Right panel: Brightness temperature changes along the jet in 2155 — 152 at 15.4 GHz. 
Filled blue circles are the measured values. Open red circles are the predicted brightness temperatures from 
shock-in-jet model with adiabatic losses dominating the radio emission assuming predominantly transverse 
magnetic field orientation to the local jet direction. The errors are given at 3(7 level. 
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4. Conclusion 

The measured sizes and brightness temperatures of the VLBI components in parsec-scale jets 
in radio galaxy 1128—047 and quasar 2155—047 are found to be consistent with emission from 
thin relativistic shocks dominated by adiabatic energy losses. The bright distinct features in these 
sources may indeed be a collection of relativistic shocks developing and expanding down the flow. 
Applying the shock-in-jet model to the VLBI observations we predict the predominantly transverse 
direction of the intrinsic magnetic field associated with the central channel of the jet in 1 128—047. 
Discrepancy between the model and measured brightness temperatures for the bright jet compo- 
nent, located at ~3.4 mas from the core in 2155—152, can be explained either by reacceleration of 
the particles due to the interaction with the ambient medium or by changing the viewing angle that 
becomes nearly zero. 

References 

[1] S.A. Baum, CP. O'Dea, G. Giovaimini, et al., ApJ 483, 178 (1997). 
[2] E.R. Craine, P.A. Strittmatter, S. Tapia, et al., ApL 17, 123 (1976). 
[3] K.I. Kellermaim, M.L. Lister, D.C. Homan, et al., ApJ 609, 539 (2004). 
[4] Y.Y. Kovalev, K.I. Kellermaim, M.L. Lister, et al., AJ 130, 2473 (2005). 

[5] A. P. Marscher, Interpretation of Compact Jet Observations, in proceedings of Parsec-Scale Radio 
Jets, ed. J. A. Zensus & T.J. Pearson (Cambridge: Cambridge Univ. Press 1990), 236. 

[6] A. P. Marscher, Variability of the Non-thermal Emission in the Jets ofBlazars, in proceedings of 
Blazar Continuum Variability, ASP Conference Series 110, ed. by H.R. Miller, J.R. Webb, and 
J.C. Noble, pp. 248-261 (1996). 

[7] R. Laing, Brightness and Polarization Structure of Decelerating Relativistic Jets in proceedings of 
Energy Transport in Radio Galaxies and Quasars, ASP Conference Series 100, pp. 241-252 (1996). 

[8] M.L. Lister et al., in prep (2009). 

[9] A. P. Lobanov, E. Carrara, J. A. Zensus, Vistas in Astronomy 41, 253 (1997). 
[10] A.P Lobanov, TP Kriciibaum, D.A. Graham, et al., A&A 364, 391 (2000). 

[11] A.P. Lobanov, J.A. Zensus, Active Galactic Nuclei at the Crossroads of Astrophysics, in proceedings 
of Exploring the Cosmic Frontier: Astrophysical Instruments for the 2F' Century, ESO Astrophysical 
Symp. Series, ed. by A.P. Lobanov, J.A. Zensus, C. Cesarsky, P.J. Diamond, pp. 147-162 (2006). 

[12] A.B. Pushkarev, D.C. Gabuzda, Yu.N. Vetukhnovskaya, and V.E. Yakimov, MNRAS 356, 859 (2005). 

[13] E. Ros, J.A. Zensus, A.P Lobanov, A&A 354, 55 (2000). 

[14] T. Savolainen, Y.Y Kovalev, A&A 489, 33 (2008). 

[15] M.C. Shepherd, Difmap: an Interactive Program for Synthesis Imaging, in proceedings of 

Astronomical Data Analysis Software and Systems VI, ASP Conference Series 125, ed. by Gareth 
Hunt and H.E. Payne, pp. 77-84 (1997). 

[16] M.P Veron-Cetty, P Veron, A&A 455, 773 (2006). 

[17] S.J. Wagner, A. Witzel, ARA&A 33, 163 (1995). 

[18] J.-H. Woo, M.C. Urry, ApJ 579, 530 (2002). 



6 



